BOREL SUMMATION OF ADIABATIC INVARIANTS 



O. COSTIN, L. DUPAIGNE, AND M. D. KRUSKAL 



Abstract. Borel summation techniques are developed to obtain exact invari- 
ants from formal adiabatic invariants (given as divergent series in a small pa- 
rameter) for a class of differential equations, under assumptions of analyticity 
of the coefficients; the method relies on the study of associated partial differen- 
tial equations in the complex plane. The type and location of the singularities 
of these associated functions, important in determining exponentially small 
corrections to formal invariants are also briefly discussed. 



1. Introduction and main results 

For many ordinary differential equations depending on a small parameter (say 
e) it is possible to construct adiabatic invariants: these are generically divergent 
expansions in e, formally constant with respect to the dynamics. Under relatively 
mild assumptions, locally there exist actual functions which are invariant to all 
orders in e and, under further conditions, even within exponentially small errors 
(see for instance ^H] and also the literature cited there). 

In the present paper we show that under suitable analyticity conditions, adiabatic 
invariants are Borel summable and their Borel sums are exact invariants in the 
regions of regularity. 

The technique that we use applies to a wide class of differential systems but for 
the sake of clarity we prefer to focus on relatively simple equations and discuss later 
how the results and methods extend. Our prototypical equation is 



(1.1) ^' + (e- 2 -V(x,e))4> = (e ->()), 

a special case of which is the one dimensional Schrddinger equation in the large 
energy limit with analyticity and decay conditions in some strip in C. 

A number of different equations in which formal invariants arise can be easily 
brought precisely to the form (|1.1|) . 

The parametrically perturbed pendulum 

(1.2) x + uj{te) 2 x = 0. (e -*■()) 
after the substitution x(t) = iv~ 1 / 2 (et)f(f^ t uj(s)ds) becomes 

Our results do not apply as such to equations with periodic coefficients but, by 
transformations, some of these equations can be brought to our setting. 
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For example, Mathieu's equation / |13| p. 404)~ in a singular perturbation regime 

e 2 h - (acos(2x) + b)h = 
can be transformed by taking h = sin -1 ' 2 (2a;) f(a cos(2a;)), a = acos(2a;), to 

2 1 / 2 a( r + ba 2 2a 2 + a 2 \ 
5 A\a{a 2 -a 2 ) (a 2 - a 2 ) 2 J J 
Also, as it will become transparent, the equations for which the methods in this 
paper apply can be higher order and/or contain nonlinear terms. 

We let 

(1.3) A = e -1 . 

Eq. (|1.1|) admits the formal solutions 1 

OO OO 

A + (\y x * Y / g k>+ {x)X- h + A_(\)e- iX * Y,9k,-W k 

k=0 k=Q 

with 



5o, ± = 1, gi± = TTji J V(s)ds, 



We consider (|l.lfl on a bounded open interval I C R with initial condition 
prescribed at some xq G I. Without loss of generality we may take I = ( — 1, 1) and 
x = 0. 

We also assume that 

#)) = /(A), 

where / is the sum of a convergent or Borel-summable series. 

Throughout this paper, Borcl summation is understood in the following way: 

Definition 1. A Borel-summable series y := yk^~ k , K G Z is a formal 

power series with the following properties 

(i) the truncated Borel transform Y — By := X)fe>o (fc^TT ^ °f V ^ as a 
nonzero radius of convergence, 

(ii) Y can be analytically continued along [0, +oo) and 

(hi) the analytic continuation Y grows at most exponentially along [0, +oo) and 
is therefore Laplace transformable along [0, +oc). 
The Borel sum y of y is then given by 

o 

(1.4) y = CBy:=Y J V^ k + ^ 

k=K 

where the sum is understood to be zero if K > and C denotes the usual Laplace 
transform. 

Remark 1. We note that although we require conditions in a C— neighborhood Af 
of [0, oo) rather than in a sector, uniqueness of y follows since Y is analytic in M 
and is uniquely defined near zero by By. 

Our main results are the following : 



^This type of expansion is best suited for Borel summability in our setting. 



BOREL SUMMATION OF ADIABATIC INVARIANTS 



3 



Theorem 2. Assume that V is analytic in e and x: 

oo 

(1-5) V(x,e)=J2Vk(x)e k 

with Vk analytic in the strip S = { — 1 < < 1}, real valued on the real line and, 
for some B, K,S > 0, satisfying 

(1.6) sup (1 + \x\ 1+s )B- k \V k (x)\ = K < oo. 

k<EN,xeS 

Then given A = e~ l > large enough, the general solution of fl.lp in (— 1, 1) can 
be written in the form 

$ = de^^+ix- A) + C 2 e~ lXx ^ (x; A) 

where </>+ and 0_ are conjugate expressions of each other and are the Borel sums of 
their asymptotic power series (J2ak + (x)X~ k and ^afc_(ir)A~ fc respectively). The 
initial condition is taken to be 

(1.7) V(0) = /(A) = 1/A. 

Furthermore, (f>+ (and (/)-) are uniquely determined in the following sense : if for 
some neighbourhood V of and some v > 0, 

<f>(x;X)=£( X (x,-)), 

where X (-,t) £ C 2 (V) for t G (0, oo) and X (x,-) £ L x v := L 1 {R+, E - Ut dt) for x G 
(—1,1) , and ife lXx (f> solves if-?. If) , on some neighbourhood of then 

(j) = tf> + . 

Theorem 3. Consider a V as in Theorem and take ip to be a solution of 
f?~^ ) in a neighbourhood ofO, which can therefore be written as 

tf, = C ie lXx dp+{x- A) + C 2 e- tXx cj)^{x- A). 

Then, in the region where the assumptions are satisfied, C(x; A; tj), ij)') :— C\C 2 is 
an exact invariant of and is the Borel sum of an adiabatic invariant 

(1.8) C = Y,Ck(x\^,ip')\- k , i-e. 

C(x;X;tp,ip') = CBC. 

For example, in lfl.2fl . in the regions where the analyticity assumptions are ful- 
filled there exists an actual invariant of the form C(t; e; x, x). C is analytic in a 
sector in e > 0, analyzable in the given sector at e = (in this case, it simply means 
that the Taylor series of the function is Borel summable to the function itself), and 
it is straightforward to check that to leading order it assumes the familiar form 

C~E/{2u 2 )- E=±(x 2 + cj 2 x 2 ) 

Remark 2. All of our results still hold if instead of (11.71) . one assumes only that 
"0(0) = /(A), where f is the sum of a Borel summable series. This is a simple 
consequence of the fact that equation (|l.lfl is linear and that Borel summable series 
form an algebra (see Lemma\^). 
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Remark 3. All of our results also hold if instead of (|1.6|) . the coefficients Vk satisfy 
the following weaker growth assumption at infinity: 

\V k {x)\ < B k g{\x\) for some B > and g e L X (R). 

2. Proof of Theorem El 

The core of the proof is contained in the case where V is independent of e and we 
start with this case. In the following C denotes any constant the value of which is 
not significant to the analysis. 

2.1. Formal derivation of a fixed-point equation. First, we (formally) manip- 
ulate (|1 . 1|) to obtain a fixed point problem. To leading order, the solutions of 
are e lXx , e~ tXx so that we look for solutions of the form 

(2.9) V = ^+ = e iXx cj> + (x) and ifi = V- = e -iAa: <£_ (x). 
Then </>+ solves 

(2.10) </>" + 2i\(j)' - V(x)(p = 0. 

We then seek for solutions of (|2.1U|) of the form (f>(x; A) := C(\(x, •)), where C is 
the usual Laplace transform. One can easily check that 

^ = C(V(x, •)) where 9(x, t) = / xfo r) dr, 

so that dividing (|2.10|) by A leads to 

f V xx + 2i* ai =V(a:)tf 



o 



(2.11) 



*(x,0) =0 

*(0,t) = / dr = i 

JO 



To control this equation, we pass to bicharacteristic coordinates, namely we let 
(2.12) s = -2ix + t 

and obtain for Q(s,t) := ^(x,t), 

U st = \v{i{s-t)/2)$> 
( 2 - 13 ) { 9(s,0) =0 

I *(M) =t 

Integrating 12.13fl yields 

(2.14) *(«,*) = J(#)(«,t), 

where 

(2.15) 



J(#)(M) =*-7 / / ^(*(si -ti)/2)$(s 1 ,t 1 )dsi* 1 



4 Jo 7t 

t + t(s-t) I I V[i((l-a)s + (a-f3)t)/2]$[(l-a)s + at,f3t] dad(3. 
Jo Jo 
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2.2. Solving the fixed-point equation: analyticity and exponential bounds. 

It is useful to note the relation between and $: 

= ®(-2ix + t,t) or $(s,i) = *(i(s -t)/2,t). 

For ^ > and f2 an open set of C 2 , let now 25^(11) be the set of analytic functions 
F over fl with finite norm 

(2.16) \\F\\„ = S up\F(s,t)\e-»W. 

n 

Let 

P ={(«,*) : |s|<P,|t|<P} 

and 

Si ={(«,*) :»(«),»(*) >0, max(|3f(t)|,|9(js)|,|9(t-a)|) < 2} U P. 

Lemma 4. Assume V(x) is analytic in S = {\x\ < R} U {a; : |3?(x)| < 1} and 
\V(x)\ < K(l + Ixl)- 1 - 5 in S for some 8 > 0. Then 

J : ?B v (Si) — » 95i/(5i) 

is a contraction for v > Zarge enough. 

Proof. First we check that J is well-defined on 25 „ (Si). Clearly, 

(s, ()eP=^ i((l - a)s + (a - /?)*)/2 G {x : |x| < P} 

for a,/3 G [0,1]. Also 

|«(*((1 - a)s + (a - /3)t)/2)| = |(1 - a)3f(s) + (a - /3)9(i)|/2 < 1, 

provided this inequality holds for (a, /3) = (0, 0), (0, 1), (1,0), (1, 1) since the above 
expression is linear in (a, 0). The condition is thus 

max(|9f(t)|, |3(s)|, |9f(t - s)|) < 2, i.e. (s,i) G S x . 

It follows easily that J maps analytic-over- Si functions to themselves. Also, \\t\\ v < 
v~ x . We are left with proving that J is contractive (in a nonlinear case, one would 
restrict J to a ball of radius Cv^ 1 ) and this follows from the following estimate 



< 



V{i{s 1 -t 1 )/2)F(s 1) t 1 )ds 1 

C\\F\\ v e v ^\ I (l + \u\)- 1 - 5 du = C\\F\\ v e^, 



(2.17) 



whence by integration 

□ 

Applying Lemma^with R = 1 and Picard's fixed-point theorem, we thus obtain 
a solution $ G 33* (St) of QFTQ . 
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2.3. Solving the original equation. 

Proof. This part of the proof is standard but we include it for convenience. To 
go back to the original equation, we just need to reverse and justify the transfor- 
mations that lead us from to H2.14fl : since $ £ fB v (Si), \& is analytic in a 
neighbourhood of (—1, 1) x (0, oo) and 

|*(M)| <Ce v ^. 

Clearly, 'J satisfies ()2.11(l . We claim that given any xo £ (—1, 1) there exists r > 
such that each partial derivative of '5 is exponentially bounded on (xq — r,xo + r) x 
(0, oo) so that (j> = X£.( i S(x, ■)) is well-defined and will solve l|2.10|l with the initial 
condition <^>(0; A) = 1/A. Indeed, 

V x = -2i$ s , ty xx = -4$ ss 

* t = $ s + $ t , * ai = -2i(<P ss + $ at ) 

But since $ is analytic in U = P U Si, letting so = — 2ix + to for x £ (xq — r, xq + r) 
and to £ R + , we have that (so,io) S U and by Cauchy's formula, 



|$ a (ao,*o)| = (2tt) 



$(s,t )/(s - s ) 2 ds 

S(s ,r) 



< C r e 



"l*ol 



where S(sa,r) denotes the circle centered at so of radius r = (1 — |xo|)/2. The 
bounds on the other derivatives can be obtained similarily and we therefore take 
Laplace transforms in H2.11(l to obtain a solution of p.l|l . l|1.7|l on (—1,1), of the 
form 

il>(x;\) =e iXx (j> + = \e iXx £(t>(x r )) where \V(x, t)\ < Ce u ^ . 

Working with 0_ = <p + (as defined in (|2.9|0 . we obtain similarily a solution of the 
form ijj{x]X) — e~ lXx 4>-. 

In the following we will rely on the fact that Borel summable series are closed 
under algebraic operations. 

Lemma 5. Borel summable series form a field. 

Though rather straightforward (see 0, ffHj; [121 )■ we provide a proof for 

convenience of the reader. The fact that Borel summable series form an algebra is 
shown in p^. Given a series y, we want to construct its multiplicative inverse (y) -1 - 
Up to factoring out a monomial dk^ K1 in the expension of y, we may always 
assume that 

y = A(l + /) for some Borel summable / = o(l). 
The inverse g of y must then satisfy 

(2.18) (1 + f)~g = 1/A. 

If F = Bf is the Borel transform of /, defined on a fixed (star-shaped-about-the- 
origin) neighbourhood £1 of [0, oo) and if G denotes the Borel transform of the 
inverse g we are looking for, we must have 

(2.19) G=1-F*G, 
where 



f F(q)G(t-q)dq:=t f F( 
Jo Jo 



(2.20) (F*G){t)= F(q)G(t- q)dq:=t F(at)G({l - a)t)da 
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Given A > and v > 0, consider the norm 

(2.21) ||F|| :=sup(A+|t|) 2 |F(t)| e -^, 

te.fi 

denned for F 6 B Vt A(Q), the space of (exponentially bounded) analytic functions 
over Q, equipped with the above norm. H2.19|l will have a solution G provided 
G — > F * G is a contraction in that space. Now, 

r\*\ 

(2.22) \F*G\<\\F\\\\G\\e v \ t \ I {A + u)~ 2 (A + \t\ - u)~ 2 du 

Jo 

= \\F\\ ||G||e^l ^"' /2 + ^ (A + U )- 2 (A + |t| - u)- 2 du 

and 

/•|*l/2 roo n 

\ (A + u)- 2 {A + \t\ - u)- 2 du <{A+ |*|/2) 2 / (A + u)- 2 du <--(A+ \t\y 2 . 
Jo Jo A 

Working similarly with the last integral in (|2.22(1 . we obtain that 
(2-23) ||F*G||<^||F||||G||. 

Hence, fixing A > large enough it follows that G —> F * G is contractive. 

Taking Laplace transforms and applying Watson's lemma [Hj, we obtain the 
(unique) solution g of H2.18JI defined by 

£G~g. 

□ 

2.4. Uniqueness of (f> + . Take <fi = £(x(x, •)) and V a neighbourhood of the origin 
as in the statement of Theorem As in Step 2.1, using the same notations to go 
from (p to $, it follows that <E> solves l|2.14|l . so that we only need to prove that J 
is a contraction for the norm 



/■oo 

F\\ Y = sup / \F{x,t)\e- vt dt 

x Jo 



and observe that since L\ C L^, when v < u' , we may choose v > as large as we 
please. Rewriting H2.15fl in terms of \ £ Y, we have 

(2.24) K{*)(a,t) :=-4(J($)(M)-i) 

= t(s-t)l I J V(x a/3 )x(x a i3,T)dT da d/3, 

where 

i 



o Jo Jo 



So 



(2.25) \K{$)\<\\ X \\Y\t{s-t)\ / / |y(a;^)|e^ t dad/3<C I /- 1 ||x||y, 

Jo Jo 

where we used <|2.17ll in the last inequality. 
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2.5. -0+ an d ip- are linearily independent. Observe that ip+ and ip- defined in 
(12.9(1 are conjugate expressions of each other and arguing by contradiction, suppose 
that for some A = A(X) £ S 1 

-0+ = Alp- 
Let B(x; A) := </> + /0_. Then 

(2.26) B = Ae- 2lXx . 

By the Lemma El B defined by 1(2. 26|) is asymptotic to a Borel summable power 
series and since <^+(0; A) = 0_(O; A) = 1/A, we have for x = 0, 

(2.27) A(X) = B(0; A) = 1 + o(l) as A -> oo, 

whereas for any other fixed value of x ^ 0, there exist if = K(x) S Z and = 
6(x) G C such that 

A(A)e- 2jAa: = B = A A '(6 + o(l)). 
Combining this equation with ((2.27(1 . we obtain 

{1 + o{l))e- 2lXx = X K (b + o(lj), 

which is clearly impossible. 

2.6. The case (|1.5H . We explain here how to adapt the argument for a potential 
V = V(x,e) depending on t = 1/A. According to 1(1.5(1 . if <j){x) = C(x(x,-)), the 
inverse Laplace transform of V<f> is given by: 

(2.28) c-W^^Voix^ + f^^^-t^^x, 

fc=i ^ >' 

where * is the convolution product defined in 1(2.20(1 . Hence, taking as before 
^(Xjt) = J x{x,t) dr, we obtain 

Following the analysis of the e-independent case, we end up with the following 
operator : 

j(*xM) = mi*,*) - zE /* r ^ti^'^ ^-^^uh) d S1 a u 

k>0 Jo Jt { )■ 

where J is defined bv ((2.15|) . Again, we must prove that J is contractive. To do so, 
instead of ((2.16(1 . we use the norm || ■ || defined in ((2.21(1 . For A large enough, it 
follows from (|!T23|) that 

(2.29) ||F * G\\ < \\F\\\\G\\ holds for all F, G analytic in ft. 

In order to control || J||, we perform a few side computations. We first estimate 

1| < C sup {t k + x e- vt + t k - x e- vt ) = C 

teR+ 

By Stirling's formula, it follows that 

(2.30) Hf*- 1 !! < C(k- 1)! k 3/2 u- k+1 . 



m — + ^ ■- 
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Next, we estimate for t > the following quantity : 

On the one hand 

while on the other hand 

f e Hti-t) ( All) dh < C f e v ^~^ dh < Civ. 
Jt/2 \A + tJ J t/2 

Hence, 

(2.31) e~ ut {A + tf f e^CA + ti) -3 *! <C/v. 

Jo 

We can now estimate || J($)||. Instead of (|2.17(l . we have 

f-t />00 

\ V k (i(s 1 -h)/2)[t k - 1 *^}( Sl ,t 1 )ds 1 <C\\^- 1 ^\\e v ^\A+\t 1 \)- 2 B k (l+M)- 1 "^ 

J s J —oo 

<Cv(-\ (k-l)l &/ 2 e v M(A+\ti\)- 2 \\§\\, 



where we used (|2.29[) and i|2.3U[) in the last inequality. Whence by integration and 
by (|2.31() . we finally obtain 



fc>o 



\\jm\<c[^ i + J2k 3/2 (^) k ) \\n<-m- 

□ 



VI V 



2.7. Further generalisations. It is not difficult to modify the proof given to allow 
for higher order equations, which, possibly after transformations, have relatively 
simple bicharacteristics, or for nonlinear dependence on ip which does not affect 
the highest derivative. In the nonlinear case, the initial condition must obviously 
be left in the form of a general sum of a convergent or Borel summable series. The 
strategies of dealing with nonlinearities are described in [7] and UJ. 

3. Proof of Theorem 01 

This follows easily from the explicit expression of C, Theorem El and Lemma 
Indeed, differentiating ip (with respect to x), we obtain 

</, = de tXx <j) + + C 2 e- lXx 4>„ 
4)' = C x e iXx A + C 2 e- lXx B 

where 

A = cf), + i\<j)+ and B — <j>_ — i\4>- 
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So that 



4 


<t>- 




4>+ i> 




B 




A 










2 




<f>+ < 


b- 






A 


B 





C := dC 2 = 



By construction, 0+ ~ ^2ak(x)X k and is uniquely determined so that using 
Lemma |3 we may conclude that 

where all considered power series are Borel summable. □ 

4. Discussion of singularities of <f 

So far we have assumed that V was analytic. If this is not so, singularities of W 
can originate in the singular points of V. We restrict the analysis to the relatively 
common situation where V has a branch-point at some point xo (we may assume 
without loss of generality xq — 0) of the form 

V(x) = x~ Vi(x), 

with V\ analytic at 0. We assume < < 1 (for different 5R/3, the analysis can 
be done similarly). 

We let P = {(s,t) : \s\ < 4^, \t\ < r} and consider the region 

Si = {(s,t) : 5ft(s) > &(t) > 0, max(|Q?(i)|, |Qf(s)|, - s)\) < 2} n P 

Assume that (s,t) — > Vi(i(s — t)/2)) is analytic in Si and that V satisfies the 
following estimate (similar to that in Lemma0J) : for some S > 0, K > 0, 

\V(x)\ < Kil + lxl)- 1 -* if |x|>l. 

Proposition 6. There exists a solution $ G C (Si) of (|2.15|) . Furthermore, letting 
vo — Vi(0), $ satisfies 

, s f \ v r s 2 ~P - (s - t) 2 ~P (s - t) 2 -P 

(4.32) »(..') = ' + 7 {. (1 _ wa ( _ w ;_ g) + (1 '_ g)( i,_ g) ' 



+ 



(1-/3) (2-/3) (3-/3) 
as i — > 0. 



}+0(t 3 ), 



Remarks 

(1) In particular, it follows easily that for s ^ 0, $ does not extend analytically 
at t = 0. Indeed, assuming the contrary, we would have Ct 3 ^ 13 = A(t) + 0(t 3 ) with 
A analytic; as t — > this forces Ait) ~ Ct?~~P . But ^4 is analytic and we must have 
for some n € N and C" 6 R, A(i) ~ C"i", which is a contradiction. 

(2) Proposition provides information about the analytic continuation (still de- 
noted <&) of the solution to (|2.15|l . Indeed, we can work as in Lemma 0] to show 
analyticity in a region where t > r > 0. Lemma0]and Proposition[|)]provide unique- 
ness in the space of exponentially bounded analytic functions over the corresponding 
region, so that by obvious imbeddings the analytic continuation $ coincides with 
the fixed point of I2.15fl . Changing variables, i.e. going back to 'J, we conclude 
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that '5 is nonanalytic at t = 0. This implies in turn a Stokes transition on cj> (and 
thus the adiabatic constant constructed in Theorem 3). 

(3) Furthermore, the fixed point procedure that leads to relation i|4.32[l can easily 
(and rigorously) provide more detailed information about the singularity manifold; 
it suffices to construct the space of functions in which the fixed point equation 
(|4.33|l below is considered in such a way that the norms entail information about 
the singularity type to be proved: the details are quite straightforward but at the 
same time quite long and we will not elaborate on them in the present paper. 
Similar constructions can be found in 0. 

Proof. From Ij2.15|l we have 

(4.33) *(«, t)=t-\ f f YM&ZW ^ ti)d8xdt 

The fact that J is defined and contractive with norm 0(t) on the functions defined 
in Si with the sup norm, follow as in Lemma ^ 

We have (1 - J)$ = t and thus for small t, $ = t + Jt + 0{J 2 t) =t + Jt+ 0{t 3 ) 
and, with vq = V\ (0) the result follows. □ 
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